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Recent results 0, [H, [3] on the stability of the spectral gap under general perturbations for 
frustration-free Hamiltonians, have motivated the following question: Does the entanglement en- 
tropy of quantum states that are connected to states satisfying an area law along gapped Hamiltonian 
paths, also satisfy an area law? We answer this question in the affirmative, combining recent ad- 
vances in quasi-adiabatic evolution and Lieb- Robinson bounds with ideas from the proof of the ID 
area law 1121. 
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I. INTRODUCTION 



Over the past decade, there has been a rapidly growing interest in the role of entanglement as a measure of 
complexity in simulating properties of quantum many-body systems Il8j. From a theoretical point of view, the 
question of how the entropy of entanglement scales with system size, has generated some spectacular results, such 
as Hasting's area law for one-dimensional gapped systems (see also [13] for partial generalizations) , with a 
stronger bound for groundstates of frustration-free Hamiltonians given by Arad et al. in as well as the recent 
i— I- breakthrough by Brandao et al. in |3(, where it is shown that exponential decay of correlations implies an area law 
,-Ch ' for the entanglement entropy of one-dimensional states. 
^4 Here, we prove that the amount of entanglement contained in states connected along gapped, local Hamiltonian 
paths satisfies a certain stability property. Specifically, we show that groundstates of gapped Hamiltonians whose 
entanglement spectrum decays fast enough to imply an area law for the entanglement entropy, are connected via 
gapped, local Hamiltonian paths to states that satisfy a similar area-law bound. The techniques used to prove this 
result apply equally well to showing the stability of the area-law for any eigenstate of a Hamiltonian with fast enough 
decay in the entanglement spectrum, as long as the connection is along a gapped, quasi-local Hamiltonian path. As 
a consequence, groundstates of gapped, frustratio'n-free Hamiltonians that satisfy the conditions for stability of the 
spectral gap under weak perturbations, found in 0, [H, [Hj] , become central objects in the study of entanglement scaling 
in 2D and in higher dimensions. This follows from the fact that the condition of local indistinguishability necessary for 
the stability of the spectral gap, also implies a rapid decay in the entanglement spectrum of the gapped groundstates 

Ov M- 

^sO \ The main result appears in Theorem [1] which makes use of Lemma [2] to convert a bound on the decay rate of 
\Q • the entanglement spectrum into a bound for the entanglement entropy of the states we are interested in. The main 
technical tool appears in Lemma [TJ where we show that the unitary evolution of the initial groundstate along a 
s ! • gapped path can be approximated, up to rapidly-decaying error, with a product of three unitaries, two of which act 
on complementary regions and the third acts along the boundary separating the two regions. We expect that such 
. a decomposition will find a variety of further uses, as it provides a more rigorous view of adiabatic evolution as a 
quantum circuit of local unitaries with depth dictated by the accuracy of the approximation we are willing to tolerate. 



II. QUASI-ADIABATIC EVOLUTION FOR GAPPED HAMILTONIANS 

We begin by introducing the central technical tool used in bounding the spread of the entanglement as the initial 
state is adiabatically transformed. The main idea is to simulate the true adiabatic evolution with another evolution 
that satisfies two important properties: 

1. The simulated evolution, introduced in Q as quasi-adiabatic continuation and further developed in (3. ITol. [T3. fl7| ■ 
is indistinguishable from the true adiabatic evolution, if we restrict our attention to the evolution of uniformly 
gapped eigenspaces. 

2. Unlike the true adiabatic evolution, the quasi-adiabatic evolution is generated by quasi-local interactions and, 
hence, transforms local operators into quasi-local operators. 
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First, let us define precisely which families of Hamiltonians we are considering. 

Definition 1. Let H(s) = Hq + J^ueA Vv,{s), where A C Z d and set b u {r) = {v € A : d(u,v) < r} to be the ball of 
radius r, centered on site u G A. Then, the following assumptions hold: 

i. H a = Y^ueA Qu> with Q u supported on b u (ro), for r > and \\Q U \\ < J\, for u G A. 

ii. V u {s) has support on b u (ro), with V u (0) — and ||V^(s)|| < J\, \\d s V u {s)\\ < J2, for s G [0, f] and u e A. 

Hi. H(s) has spectral gap 7(3) > 7 > 0, with differ entiable groundstate \tjjo(s)) , for all s G [0, 1]. In particular, Hq 
has unique groundstate \ipo(0)} :— |?Ao); with spectral gap 7(0) > 7. 

For the above one-parameter family of Hamiltonians H(s), there exists a family of unitaries U s @, satisfying: 

\M»)) = u.\ih), d s u s = iv s u s , f/ = i, v«g[o,i], (1) 

where the generating dynamics, 2? s , has the following decay properties: 

V s = J2 ^( u ; r )' supp(2? s (u;r)) = 6 u (r), \\D s (u; r)\\ < 2 J 2 / 7 (r - r ), (2) 

uEA r>ro 

for a sub-exponentially decaying function / 7 (r) (e.g. exp{— c rj In 2 r}, for c > 0), with / 7 (0) = 1 and decay rate 
proportional to 7/wo, where vq ~ Ji r-Q -1 is the Lieb-Robinson velocity for the family of Hamiltonians {H(s)} s >o [ill, 



III. DECOMPOSING THE QUASI- ADIABATIC EVOLUTION 

The above decay estimates on the generator of the unitary U s , allow us to decompose the quasi- adiabatic evolution 
into a product of three unitaries, two of which act on disjoint subsets of the lattice and a third one coupling the 
disjoint evolutions at the boundary of the two sets. More importantly, the error in approximating the true adiabatic 
evolution decays sub-exponentially in the thickness of the boundary chosen for the coupling unitary. 

Lemma 1. Let U s denote the unitary corresponding to the quasi- adiabatic evolution of a one-parameter, gapped family 
of Hamiltonians H(s), as defined above. For A C A, set: 

1. I A (R) = {x G A : d(x, dA) < R}, 

2. E A {R) = {x G A c : d(x, dA) < R}, and 

3. 8A{R) = I A (R) U E A (R). 

Then, there exist unitaries V S (A) , U S (A C ) andU s (8A(2R)) with non-trivial support on A, A c and dA(2R), respectively, 
such that the following bound holds: 

\\U S (A) ® U S (A C ) U s (dA(2R)) - U.\\ < e s (R) := Cl (e c ^'^ - l) \3A\ / 7 (c 3 i?), (3) 

for dimensional constants ci, 02,03 > 0. 

Proof. We sketch the proof, which uses the Lieb-Robinson bounds developed in 0, [l(|. First, define the following 
unitary operator: 

V S (A) = UlU s (A)®U s (A c ), (4) 
with d s U s {A) = iV s {A)U s {A), U A {0) = 1 and d s U s {A c ) = W S {A C )U S (A C ), U A c(0) = 1, where we define, for Z C A: 

V S {Z)= J2 V s (u;r), (5) 

«GA, r>0: b u (r)CZ 

noting that supp(2? s (Z)) C Z and, hence, supp(£/ s (AT)) C X. Now, we set: 

F s (dA) = V s -{V s {A)+V s {A c ))= V s(u;r), (6) 
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where X n dA ^ <^=> (lni^0)A(ln4 c ^ 0). The operator T s {dA) can be well approximated by the locally 
supported V s (dA{R)). To see this, note that: 

F s {dA)-V s {dA{R)) = Y J E v s{u;r), (7) 

ueAr>d(u) 

where d(u) — max{i? + 1 — d(u, dA), 1 + d(u, dA)}. We may partition A as follows: 

A = U k >oB k , B k = {ueA: d(u, dA) = k}, (8) 

noting, further, that \B k \ < T,uedA(\ b u( k )\ - K(k-l)\) < c d k d ~ l \dA\, for k > 1 and |fl | = \dA\, where c d = 2d can 
be thought of as the area of the unit ball in Z d . This implies that: 

\\F s (dA)-V s (dA(R))\\<Y^ £ E WV'MW (9) 

k>0uGB k r>d{u) 

/L(-R+1)/2J \ 

<2J 2 E ^ E Mr~ro)+ E W E (10) 

\ fc>0 r>\(R+l)/2\ fc>l+L(-R+l)/2j r>l+k J 

<2J 2 \dA\((l + (R+l) d )F 1 (\(R + l)/2-]-r )+c d E k^F^l + k - r Q )) , (11) 

fc>l+L(-R+l)/2j 

where F 7 (s) := E r >* A<( r )- 

At this point, given the decay rate of / 7 @, it should be clear that for dimensional constants Co, do > 0: 

||.F a (0A)-X> a (&4(JZ))|| < Co (J 2 /7)|aA|(l + J R) d+do /7(r(^+l)/2l -r ). (12) 

By differentiating both sides of (gj), we get: d s V s (A) = -i T s (dA) U s ) V„(A), V (A) = 1. We may approximate 
the unitary V S (A) with W S (A) generated by: d s W s (A) = -i (W s V s (dA(R)) U s ) W S (A), W (A) = 1, such that: 

\\V S (A) - W S (A)\\ = \\Wt{A)V s {A) ~ 111 < T \\ d tWkA)V t {A)\\ dt < \s\ sup \\F t (dA) - V t {dA{R))\\. (13) 

Jo te[o,s] 



Combined with (| X2[) . the above bound implies: 

\\V,(A)-W a (A)\\ <Q,(J 2 / 7 )H|M|(l + i2)«^*/ 7 (r(i2 + l)/2l-ro), c ,d > 0. (14) 

Finally, we approximate Wg(A) with the unitary W s (dA(2R)) given by: 

d s W s (dA(2R)) = * (l/2(0A(2.R))X> a (0A(iZ)) E7.(0A(2i2))) W s (dA(2R)), W (dA(R)) = I, (15) 

with the unitary U s (dA(2R)) generated by V s (dA(2R)). Following (fl3"]l. we have: 

||W S (A) - W s (9A(2fl))|| < |s| sup \\u}v t {dA{R))U t - U}(dA(2R))V t {dA(R))U t (dA(2R))\\. (16) 

te[o,s] 

At this point, we cannot use the Lieb- Robinson bounds developed in [Tolj . [3. Thm. 4.6] directly, since we are dealing 
with evolutions of the form UJ Oa U s instead of U s Oa UJ. Nevertheless, by setting: 

F t (O x ) = U\ O x U t - u}{dA{2R)) O x U t (dA(2R)), (17) 
for an operator Ox with support on X C A, we get: 

F t (Q x )= f (Uji[O x ,V s -V s (dA(2R))}U s + F s (i[O x ,V s (dA(2R))]))ds, (18) 
Jo 

which implies: 

\\F t (O x )\\<\t\ E SU P \\[OxMu;r)]\\+ E sup \\[O x ,V s (u;r)]\\ f \\F s (O x (u;r))\\ds, (19) 

b u (r)CdA(2R) se IM b u (r)CdA(2R) sE [°>*1 ^° 
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where Ox{u; r) — [Ox,D s (u; r)]/|| [Ox, T> s (u; r)]\\. At this point, we may use the recursive argument found in [Tll.[l6|. 
setting Ox — T> t (dA(R)) and recalling the sub-exponential decay of ||2? s (u; r)||. Setting: 

e s (R) := c x (>(./ 2 /7)M _ 1) \ dA \ f y (c 3 R), (20) 

for dimensional constants ci,C2,C3 > 0, we get from (fT4|) and (fTG)) : 

||C/ S - E^)0.(A c )Wj(M(2J2))|| < ||y s t(A) - + - Wt{dA(2R))\\ < e s (R). (21) 

Renaming U s (dA(2R)) as Wj (&4(2.R) ) completes the proof. □ 



IV. BOUNDING THE ENTANGLEMENT ENTROPY: 

At this point, we are ready to apply Lemma [T] to study the entanglement entropy of states adiabatically connected 
to an initial state \ipo), whose Schmidt coefficients across a cut A : A c satisfy a certain rapid-decay condition. In 
particular, we assume that for a rapidly-decaying function f A (-), with f A (0) = 1 and A a convex subset of A: 

]T a (a) < f A (R), ReN, (22) 

a > N R\dA\ +1 

where N is the maximum single-site dimension and \ipo) = J2 a >i V a o( a ) IV^o^)} ® l^.ofa)) , the Schmidt decom- 
position of \ipo) across the cut A : A c , with Schmidt coefficients in decreasing order. 

The proof follows an argument similar to the one found in First, we approximate the initial state l^o) with 
the family of states {\iPo,r)}r>o, where: 

|^o,r) = -p= V* V a o(a) \ip A ,o( a )) ® \ipA°,o( a )) 

and cr — J2 a >i a o( a ) — 1 ~ /a(R), where we used (|22|) . 

The next step is to construct states with bounded Schmidt rank and increasing overlap to the adiabatically evolved 
state \4>q{s)) = U s \ipo). To accomplish this, we define the family of states: 

\ihM) = U M) ® U S {A C ) U s (dA(2R)) \^ R ) . 

Setting the overlap P(R) := | (ipo(s), 4>o,r(s)) | 2 to be, say, at least 1/2, fixes the minimum boundary thickness Ro we 
consider in our approximation. Then, we use Lemma [5] to show that the entanglement entropy of \ipo(s)} satisfies a 
non-trivial bound for all s € [0, 1]. 

Theorem 1. For the gapped family of Hamiltonians H{s) of Definition^ let \ipo(s)) denote the groundstate of H(s). 
Then, if the groundstate |^/)q(0)) satisfies the decay condition h2ty) , the entropy of p s {A) = Tr^c \ip (s))(?po(s)\ is 
bounded as follows: 

S(p s (A))<b(l + ci)R \dA\+hi, R = mm{ReN: f A (R Q ) + 2e s (R ) < 1/2}, (23) 
where c\ = J2n>i n ^( n ) an< ^ — — J2n>o^( n ) ^ n 3( n )> with 

S(n) = (f A (nR ) - f A ((n + l)Ro)) + 2(e s (nR Q ) - e s ((n + l)R )) 

and e s (-) defined in\3[ 

Note that c\ is a constant, as long as 5(n) decays faster than n~^ 2+e \ Moreover, h% = — Y^,n>i ^( n ) ^ n S(n) < c% lni, 
for 5(n) > i"", n > 1, t > 1 and hi < 1/lnt if 6{n) < t~ n for n > 1 and t > 1. 

Proof. Fix s £ [0, 1] and let p s {A) = Tr^c \ipo(s))(ipo(s)\. Since we will be using the Schmidt decomposition of the 
state \ipo(s)), let us introduce it here: 

IV>o(s)) = V^W \^A, s ( a )) ® lV>4<>(aO) > 

a>l 
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where J2a (J s(a) — 1 and {\ip Aia (a))} , {|V>A c ,s(a))} form orthonormal sets supported on A and A c , respectively. We 
order the Schmidt coefficients of \ipo(s)} in decreasing order such that a < /3 => &s( a ) > cr s {0). Moreover, we note 
that: 

p s {A) = YjVs{a)\il) A ,s{a)){i>AA<x)\ and S(p.(A)) = - ^ <r* (a) In a s (a). (24) 

ct>l a>l 

We define: 

P(R) = \(Ms),^o,r(s))\ 2 , A s (R) = U s (A)®U s (A c )U s (dA{2R))-U s (25) 

and use Lemma [TJ in order to relate P(R) with the sub-exponentially decaying error e s (R): 

P(R) = | (Vol 1 + U$A S (R) \^ R ) | 2 > | {^M,r) I 2 - 2||A s (i?)|| > f - (f A (R) + 2e s (R)), (26) 

noting that (ipo,i>o,n) = \/cr- Setting R to be the smallest R such that fA{R) + 2e s (R) < 1/2, the decay of Ja(R) 
and e s (R) implies: 

R>R =^ P(R) > 1/2. (27) 

Now, the state \ipo. r{s)) has Schmidt rank bounded by Ur ■ N 4R ' dA ', where k^ = N R \ 9A \ is the Schmidt rank of 
\iPo,r}- To see this, first note that the Schmidt rank of U s (dA(2R)) \ipA,o( a )) ® \^pA"fl{a)) is the same as that of 
U 8 (A) (g) U S {A C ) U s (dA{2R)) \ipA,o{a)) ® \4> 'a-= ,o(oi)) , along the boundary of A. It remains to see how the action of 
U s (dA(2R)) on \tpA,o( a )) ® |V'A c ,o(a)) affects the Schmidt rank. 

Since U s (dA(2R)) is an operator acting non-trivially only on sites in a subset of dA(2R), we have the following 
general decomposition: 

Di 

U s (dA(2R))= £ 1a\Ia.(2R) ® E(a,(3) <g> G(a,(3) ® 1{a<=)\e a (2R), (28) 

where the matrices G(a,/3) act on sites in Ea(2R) and the matrix units E(a,/3), which act non-trivially on Ia(2R), 
form an orthonormal basis for Dj x Dj matrices. Moreover, Di < N' Ia ^ 2R " < N 2R ' dA ', for convex A and N the 
maximum dimension among the single-site state spaces (e.g. N — 2 for a system of qubits). 
To bound the Schmidt rank of 

|^,o(7)) := U s (dA(2R)) |Va,o(7)) ® |^a«,o(7)) , 
we trace over sites in A and study the rank of the operator: 

Di 



Tr A (|V^,o(7))(V-B,o(7)l) = E 

a=l 

with 



J2 cifi,p)\F a {p,i)){F a {p,i)\ 



(29) 



\F a (j3,j)) = G(a,/3) ® 1a^)\s a (2H) 1^,0(7)) 

and 

c(/3,/3') = (Va,o(7)I Ia\i a (2r)®E{(3',(3) \^ A ,o{l)) 



coming from ([28]) . 

Clearly, as a sum of Dj matrices each with rank at most Dj, the above operator has rank bounded by Dj and, 
hence, the Schmidt rank of U s (dA(2R)) \4>a{i,5)) is bounded above by N 4R \ dA \. Now, the final step is to note that 
if \ipo,R.) has Schmidt rank ku, then there are exactly ku states ^5,0(7)), each with Schmidt rank at most N 4R \ dA \, 
and the bound for the Schmidt rank of the approximation \i}o,R,(s)) = X^=i V a o(l) U S (A) (gi U S (A C ) |?As,o(7)) 
follows. 

The next step is to relate the overlap P(R) to the Schmidt coefficients of \ipo(s)). More specifically, we will now 
show that for R > 0: 

]T a s (a)>P(R). (30) 
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To prove this, first introduce the Schmidt decomposition of \i/)o,r(s))'- 

r r 

\ifoM) =EV^M) ® \*A°(P)) , with = 1 andr < AT 5 ™, (31) 

/3=1 ,8=1 

as we have already demonstrated. For notational convenience, let Ma{<x,0) = \(^A((3),ipA,s(ct)}\ and Ma<=(o!, /3) = 
|($A-(^),-0A% s (a))l ■ Note that since {|$a-(/3))}, {|^a, s (c*))} and are orthonormal sets, 

Bessel's inequality implies that 

r 

^M A {a,P) 2 < 1 and ^M A o(a,/3) 2 < f, 

as well as 

a a,/3 

Then, an application of the triangle inequality followed by Cauchy-Schwarz gives the following upper bound for 
P(R) = I (ifo,R(s)>M a )) I 2 : 

KVo,ii(s),Vo(s)) I 2 < ^y^RV^)^(«^)M Ac (a,« 
< [Ea s (a)M A (a,/3) 2 ( £ r(/3)M A c (a, 

The last inequality follows from Schur convexity of f([p(a)]) — ^ Q a s (a) p(a) and the observation that the vector 

[1, 1, . . . , 1, 0, . . . , 0], with at most r ones, majorizes Yla ^Ta(I) /3) 2 > Ma(2, /3) 2 , ■ • • ■ 

To see that f([p(a)]) is Schur convex, note that if we set S p (a) = 52k=iP(k) an d A(a, /3) = cr s (a) — <x s (/3) then the 
condition that [p(a)] majorizes [<?(«)] (p >: q) becomes p y q <4> S p (a) > S q (a), Va. Moreover, 

/({P(«)» - /({?(«)}) = £ A ( a ' a + !) ( 5 » - 5 ») > 

a 

which is non-negative since we have arranged the a s {a) in decreasing order so that A(a, a + 1) > 0, Va. 

Now that we have demonstrated (f5U|) . we may use it in combination with (j2"T)l to show that satisfies an 

entropy bound. Using (126H27[) and (|30l) . we have for n > 1: 

£ tr.(a) < 1 - P(nJ2o) < /A(»i2o) + 2e s {nR ). (32) 

a>AT 5 "- R o 

We now maximize the entropy 

S(p.(A)) = -'£ i <T.(a)H<r.(<x)) 

a=l 

subject to the constraint (J32J . Following the notation of Lemma[2j set s n — ]\[ 5nR o\dA\ ^ n >\ anc j f(n) = Ja^uRo) + 
2e s (ni? ), n > 1 and /(0) = 1, noting that /(l) < /(0), by our choice of i? - Then, for r = N m °\ aA \ and si = 
iV 5 ^ !^!, Lemma [2] implies that: 

S(p,(A))<5(l + ci)i2o|5i4| + /ii, (33) 
where ci = 5Z n>1 n ^( n ) and = — 5Z n>0 <5(?i) hi<5(n), with 

*(n) - (/A(ni? ) - /a((»» + l)i?o)) + 2(e s (ni? ) - e s ((n + 1)R )). 

□ 
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We turn, now, to the following Lemma, which relates the decay properties of the entanglement spectrum to a bound 
on the entropy of entanglement: 

Lemma 2 (Entropy bound). For p a density matrix, let X) Q >i a ( a )\4 , ( a )}('>P( a )\ be its spectral decomposition, with 
o~(a) in decreasing order. Assume that there is an increasing sequence of integers {s n }n>o with s n < s n+ \ < r s n , for 
n > 1 and si > sq = 0, such that the following constraint holds for a strictly- decaying function f(n), with /(0) = 1: 

E cr{a)<f{n), n > 0. (34) 

a>s„+l 

Then, the entropy of p, given by S(p) = — Yl a >i fT ( a ) m a ( a )i satisfies the following bound: 

S{p) < In si + ci In r + hi, 
where c\ = J2n>i n &{ n ) an d h\ = — J2n>o ^( n ) ^ n S(n), with S(n) = f(n) — f(n + 1). 

Proof. Since the Shannon entropy is Schur-concave, the von Neumann entropy S(p) is bounded above by the Shannon 
entropy H(~jl) of any probability distribution {/i(a)} Q >i consistent with the constraint (I34[) . that is majorized by 
{a(a)} a >i (i.e. ft < &.) 

Define S(n) = f(n) — f(n + 1), n > 0. For 1 < a < s±, the constraint J2 a >i a ( a ) = 1 implies: 



E °{a) > /(0) - /(l) =^ V p(a) = 6(0) =^ p(a) = ^1 , 1 < a < s x . 
z — ' z — ' Si 

a=l a=l 1 

Similarly, for a G [s„ + 1, s n +i], n > 1, we see that we should choose ~jt to satisfy: 



s„+i 



E ^ = /(«) - /(« + 1) m(«) = J(n) ■ 

Gathering terms, we decompose the Shannon entropy H{jt) as J^^Lo Hn{~ji), where we have defined H„(~jl) = 
~J2''a=s +i ^( a ) > 0. We can bound each H n (j}) as follows, recalling that s n +i < sir" and, hence, 

S n +1 - S„ < Sir": 

ffn(7t) < ~^(n) ln<5(n) + nS(n)\nr + S(n)\n si, n > 0. (35) 
Using 2 n >o ^( n ) = 1) we S e * the bound: 

(~tt) < In si + ci In r + fti , 

where we defined ci = Xm>i n5(n) and /ij = — J3 n>0 ln<5(rt). This completes the proof of the Lemma. □ 



V. CONCLUSION 



We have shown that states satisfying an area law are connected via gapped, local Hamiltonian paths to states that 
satisfy a similar bound on their entropy of entanglement. This result holds for lattice systems in any dimension. Since 
the above result relies only on Lieb- Robinson bounds and locality properties of the quasi- adiabatic evolution, it can be 
shown that even gapped paths with long-range interactions satisfy the bound (|23[) . as long as the decay is fast enough 
(i.e. power-law decay with a dimension-dependent exponent.) We hope this work will motivate interested readers to 
work out the details of this argument, which would have important implications for the classification of phases, given 
the connection of states satisfying an area-law to Matrix Product States (MPS) in ID and Projected Entangled Pair 
States (PEPS) in 2D 0,(3. 

Acknowledgements: The author would like to thank M. B. Hastings for helpful remarks on the use of quasi-adiabatic 
evolution in the context of entanglement entropy, as well as fruitful discussions with B. Nachtergaele, N. Schuch and 
A. Gorshkov. 
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